Crossover from Fermi liquid to Wigner molecule behavior in quantum dots 
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The crossover from weak to strong correlations in parabolic 
quantum dots at zero magnetic field is studied by numerically 
exact path-integral Monte Carlo simulations for up to eight 
electrons. By the use of a multilevel blocking algorithm, the 
simulations are carried out free of the fermion sign problem. 
We obtain a universal crossover only governed by the density 
parameter r s . For r a > r c , the data are consistent with a 
Wigner molecule description, while for r a < r c , Fermi liquid 
behavior is recovered. The crossover value r c « 4 is surpris- 
ingly small. 

PACS numbers: 73.20. Dx, 71.10.Ay, 71.10.Ca 



Quantum dots can be considered as solid-state artifi- 
cial atoms with tunable properties. Confining a small 
number of electrons N in a two-dimensional electron gas 
in semiconductor heterostructures, a number of interest- 
ing effects arising from the interplay between confine- 
ment and the Coulomb interaction between the electrons 
can be observed 10,0. Since the confinement potential 
is usually quite shallow, the long-ranged Coulomb inter- 
action among the electrons plays a prominent role, and 
in contrast to conventional atoms effective single-particle 
approximations quickly become unreliable. In the low- 
density (strong-interaction) limit, r s — > oo, classical con- 
siderations suggest a Wigner crystal-like phase with elec- 
trons spatially arranged in shells fl. With quantum 
fluctuations, such a phase is best described as a Wigner 
molecule. In contrast, for high densities (weak interac- 
tions), r s — > 0, a Fermi liquid- like description is expected 
to be valid, where it is more appropriate to think of the 
behavior as resulting from the single-particle orbitals be- 
ing filled. The noninteracting limit Q is then typically 
used as a starting point for the theoretical description of 
quantum dots. 

To date, no reliable information exists for the crossover 
between these two limits. This is mainly due to a com- 
plete lack of sufficiently accurate methods that are able 
to cover the full range of r s , especially when no magnetic 
field is present. Exact diagonalization techniques are lim- 
ited to very small particle numbers and small r s , other- 
wise a large error due to the truncation of the Hilbcrt 
space arises ||. Hartree-Fock calculations become in- 
creasingly unreliable for large r s and are known to incor- 
rectly favor spin-polarized states Similarly, density 
functional calculations JtJ introduce uncontrolled approx- 
imations in the absence of exact reference data. In prin- 



ciple, the quantum Monte Carlo (QMC) method is the 
best candidate for producing reliable data for quantum 
dots. Unfortunately, the notorious fermion sign prob- 
lem makes direct QMC simulations almost impossible [|| . 
To avoid the sign problem, the fixed-node approximation 
and a related variational approach have been employed 
in Ref. ||, but the results are no longer exact. 

In this study, we adopt a radically different approach 
to fermion QMC simulations, based on the recently devel- 
oped multilevel blocking (MLB) algorithm The 
MLB algorithm is able to provide numerically exact QMC 
results free of the sign problem. In this Letter, we re- 
port large-scale simulation results for the weak-to-strong- 
correlation crossover for up to eight electrons. The nu- 
merical results at large r s are shown to agree with a 
Wigner molecule description. 

Model — We study a two-dimensional parabolic quan- 
tum dot at zero magnetic field, 



H = 




(1) 



where the positions (momenta) of the electrons are de- 
noted by Xj (pj). The effective mass is m*, and the 
dielectric constant is k. The MLB calculations are car- 
ried out at fixed N and fixed S = (JVf - N±)/2, the 
z-component of the total spin. We present results for 
the energy, E = (H), the radial charge and spin den- 
sities p(r) and s z (r) normalized to J Q dr2irr p(r) = N 
and J dr lixr s z (r) = S, and the two-particle correlation 
function 



9s( x ) 



271-Z 2 , 




(2) 



g s is isotropic, and with y — t/Iq prefactors are cho- 
sen such that / °° dyyg s (y) = 1. The length scale 

lo = y/h/m*u)o from the confinement allows the interac- 
tion strength to be parametrized by A = lo/a = e 2 / kujoIq, 
where a is the effective Bohr radius. For any given N and 
A, the dimensionless density parameter r s can then be ob- 
tained from the data, r s = r*/a, where r* corresponds 
to the first maximum in g s (r). The values for r s 
obtained this way agree well with the predictions of an 
electrostatic point-charge model |l^] . In all simulations, 
the temperature was set to T = 0.1 Hcoo/kB- 

Method — The simulation method is based on a stan- 
dard discretized path-integral representation of the ob- 
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servables of interest, where the sampling is done accord- 
ing to the MLB algorithm described in detail in Ref . [^0) . 
A sample number of at most K = 600 was sufficient 
to eliminate bias from the algorithm, and at the same 
time cured the sign problem. The simulations have been 
carried out on up to five levels in the MLB scheme. Data 
were collected from several 10 4 samples for each parame- 
ter set {N, S, A}, with a typical CPU time requirement of 
a few days (for each set) on a SGI Octane workstation. 
As a validation for this procedure, we have accurately 
reproduced the exact diagonalization results for N = 2 
electrons ||. 
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Wigner molecule behavior. The classical shell filling se- 
quence has been predicted recently ||. For TV < 6, there 
is only one shell, but the sixth electron enters a new in- 
ner shell (1-5). Furthermore, for N — 7 the shell filling is 
1-6, and for N — 8 it is 1-7. These predictions are in ac- 
cordance with our data. Additional simulations for up to 
12 electrons at A = 8 [not shown here] further verify the 
classical filling sequence. The only exception is N = 10, 
where we find a 3-7 instead of the predicted 2-8 structure. 
Clear indications of a spatial shell structure at N > 6 can 
be observed even for A = 4, albeit significantly less pro- 
nounced than for A = 8. For A _> 4, the charge densities 
are found to be quite insensitive to S. This is expected 
for a Wigner crystal where particle statistics and spin in- 
fluence energies or density correlations only weakly. Our 
numerical results for the spin density in this regime sim- 
ply follow the corresponding charge density according to 
s z (r) ~ (S/N)p(r). A significant S'-dependence of charge 
and spin densities is observed only for weak correlations. 
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FIG. 2. Numerical results for £iv(r a ). Statistical errors 
are of the order of the symbol size. The dotted curve, given 
by exp(— r s /r c ) with r c — 4, is a guide to the eye only. The 
inset shows the same data on a semi-logarithmic scale as a 
function of *Jrl. The dashed line is given by Eq. (^) . 



FIG. 1. Density p(r) of the spin-polarized state (S = N/2) 
for A = 2 and A = 8. Dashed, solid, dashed-dotted, and 
dotted curves correspond to TV = 5, 6, 7, and 8, respectively. 
Units are such that lo = 1. 

Charge and spin densities — Figure [l] shows the charge 
density p(r) of the spin-polarized state for N = 5 to 8 
electrons. While for A = 2 increasing N does not change 
p(r) qualitatively, the situation is different for strong in- 
teractions (A = 8), signifying the onset of shell forma- 
tion in real space. Such a structure is clear evidence for 



Crossover To quantitatively investigate the 

crossover from weak to strong correlations, we employ 
the quantity 



oc 
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normalized such that in the absence of interactions £jv = 
1. The correlation function g s (r) in Eq. (^|) is a very sen- 
sitive measure of Fermi statistics, in particular revealing 
the spin-dependent correlation hole. Since interactions 
tend to destroy the Fermi surface, the spin sensitivity 
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of gg(r) is largest for a Fermi gas, r s = 0. In fact, for 
r s — > oo, the correlation function gs(r) becomes com- 
pletely spin- independent. Hence the quantity £at(Y s ) de- 
cays from unity at r s = down to zero as r s — > oo. The 
functional form of this decay is indicative of the crossover 
under consideration. 

As seen in Figure ^, the crossover curve £jv(?" s ) be- 
comes remarkably universal and depends only weakly on 
N. Its decay defines a crossover scale r c , where a simple 
exponential fit for small r s yields r c rs 4. For r s > r c , 
the functional form of £(r s ) is better described by 

£(r s ) oc exp (-y/rJrOj , (4) 

where r' c rs 1.2. We mention in passing that the WKB 
estimate for {i/is\H\i/is) — {^S^H^S') exhibits the same 
large-r s dependence. One can therefore argue that the 
spin sensitivity of the square of the eigenfunctions |i/'s| 2 
also shows this behavior, and thereby rationalize Eq. 
The value r c rs 4 is consistent with the appearance of 
spatial shell structures in the density profile. In addi- 
tion, the energy spectrum is in accordance with a Wigner 
molecule description for r s > r c (see below). Summa- 
rizing, the crossover from weak to strong correlations is 
characterized by the rather small value r c rs 4. 

Parenthetically, we contrast this result with the values 
r™ c rs 37 found for clean jl3| and rf s rs 7.5 for disor- 
dered 0] unbounded systems. The latter result provides 
evidence that breaking the continuous translation invari- 
ance stabilizes the crystallized phase. The even smaller r c 
found here should then be due to the confinement. Cau- 
tion is adviced with the thermodynamic limit, luq — > 
with r s fixed, where plasmon modes eventually govern 
the low-energy physics. There the value r™ c becomes 
relevant. For GaAs-based structures with r s = 4, we 
estimate that spin-sensitive properties loose their signif- 
icance only for very large electron numbers, N >, 10 . 

Energies — MLB results for the energy at different pa- 
rameter sets {N, S, A} are listed in Table 0. For given N 
and A, if the ground state is (partially) spin-polarized 
with spin 5, the simulations should yield the same ener- 
gies for all 5' < 5. Within the accuracy of the calcula- 
tion, this consistency check is indeed fulfilled. 

Detailed data are given in Table B. For N — 3 elec- 
trons, as r s is increased, a transition occurs from the 
5 = 1/2 to a spin-polarized 5 — 3/2 ground state at 
an interaction strength A rs 5 corresponding to r s rs 8. 
For N — 4, we encounter a Hund's rule case. By using 
perturbation theory in r s , one may show that the inter- 
actions lead to a 5 = 1 ground state. From our data, this 
standard Hund's rule (which applies for small r s ) is seen 
to hold throughout the full range of r s , and the ground 
state spin stays 5 = 1 even for large r s . A similar situ- 
ation arises for N = 5 electrons, where the ground state 
is characterized by S = 1/2 for all r a . Turning to N = 6, 
while one has filled orbitals and hence a zero-spin ground 



state for weak correlations, our results for A = 8 reveal a 
transition to a 5 = 1 ground state as r s is increased. A 
similar transition from a 5 = 1/2 ground state for weak 
correlations to a partially spin-polarized 5 = 5/2 ground 
state is found for N = 7 electrons. Finally, for N = 8, 
as expected from Hund's rule, a 5 = 1 ground state is 
observed for small r s . However, for A > 4, corresponding 
to r s > 10, the ground state is seen to have spin 5 = 2, 
in contrast to the conventional Hund's rule prediction. 

For strong correlations, r s > r c , the energy levels and 
their spin splittings differ considerably from what is ex- 
pected from a Fermi liquid-like orbital picture. In partic- 
ular, the ground-state spin 5 can change and the excita- 
tion energy of higher-spin states becomes much smaller 
than Hu!q. This level structure reflects the onset of shell 
formation in real space. In fact, our large-r s data in 
Table | can be rationalized by starting from crystallized 
electrons located at positions fixed by electrostatics, and 
then evaluating the quantum corrections due to particle 
exchange processes, namely rotation and tunneling. Gen- 
eralizing the method of Ref. |lfj by using semiclassical 
estimates and group theory to satisfy the Pauli principle, 
detailed predictions for low-energy spectra can be made 
in terms of such a Wigner molecule. While a detailed dis- 
cussion of the Wigner molecule will be given elsewhere, 
it is already apparent from our discussion above that the 
value for r s where the ground state spin changes is not 
given by r c R! 4 but is typically larger. Therefore such 
transitions should be amenable to the Wigner molecule 
concept, which is indeed the case. 

The above findings for the energy imply several novel 
and nontrivial consequences for transport experiments 
made by weakly coupling the dot to electrodes. In par- 
ticular, the addition energies following from Table || de- 
termine the positions of the conductance peaks measured 
experimentally by capacitance spectroscopy or by lin- 
ear transport P|. Furthermore, high spins occurring close 
to the ground-state energy can lead to negative differen- 
tial conductances in transport measurements, or even to 
the disappearance of a conductance peak at low temper- 
atures (spin blockade) [|l6| . According to Table H, N = 6 
would be a possible candidate in which to find negative 
differential conductances for the transition to N — 5 for 
r s > 15, corresponding to Hujq < 0.4 meV in a GaAs-based 
quantum dot. Furthermore, the spin-polarized ground 
state for N = 3 at A > 5 implies that the direct transition 
into the 5 = ground state for N = 2 is spin-forbidden. 
The corresponding conductance peak should then disap- 
pear for hujQ < 0.5 meV. A similar situation arises for the 
N = 7 to N = 6 transition at %wq < 0.4 meV. Such phe- 
nomena cannot occur in the weakly interacting regime 
r s < r c , where entering or escaping electrons are accom- 
modated in effective single-particle orbitals together with 
their spins. We note that sufficiently large quantum dots 
allowing for experimental studies of the Wigner molecule 
phase are within reach of current technology |17[]. 
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To conclude, we have presented numerically exact 
QMC results for parabolic quantum dots covering the 
full crossover from weak (r s — * 0) to strong (r s — > oo) 
correlations. The turnover from Fermi liquid to Wigner 
molecule like behavior is basically independent of the 
particle number and characterized by an astonishingly 
small crossover scale, r c ps 4. Energy spectra in the low- 
density regime r s > r c differ from single-particle expec- 
tations but can be described within a Wigner molecule 
approach. Detailed predictions have been made for this 
Wigner molecule phase, which should be directly acces- 
sible to current experiments. It is straightforward (and 
left to future MLB studies) to study other confinements 
or interaction potentials, or to include a magnetic field. 
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TABLE I. MLB data for the energy for various {N, S, A} 
parameter sets. Bracketed numbers denote statistical errors. 
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